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Abstract 

According to the widely accepted notion, the Schrodinger equation 
(SE) is not derivable in principle. Contrary to this belief, we present 
here a straightforward derivation of SE. It is based on only two fun- 
damentals of mechanics: the classical Hamilton- Jacobi equation(HJE) 
and the Planck postulate about the discrete transfer of energy at micro- 
scales. Our approach is drastically different from the other published 
derivations of SE which either employ an ad hoc underlying assump- 
tion about the probabilistic or the statistical nature of the micro-scale 
phenomena, or rely on the prior knowledge of SE and arrive at it by 
introducing a new postulate - neither present in classical mechanics 
nor following from experiments - with a suitable but physically unjus- 
tifiable choice of a key arbitrary constant. 
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" ...with regard to the origin of the presented recipe 
the Schrodinger equation we have to unequivocally 
say: there is no derivation and there cannot be 
any. To derive means to obtain it logically from 
other laws. The Schrodinger equation has been 
guessed." [1] 



1 Introduction. 

Thirty years after the above categorical statement by L.I.Mandelshtam, R. 
Feynman wrote [2], "For a particle moving freely in space with no forces, no 
disturbances, the correct law of physics is 



Where did we get that from? Nowhere. It is not possible to derive from 
anything you know. It came out of the mind of Schrodinger, invented in 
his struggle to find an understanding of the experimental observation of the 
real world" . 

In turn, the "inventor" himself, E. Schrodinger, admitted in the second 
of his seminal papers [3]: "So far we have only briefly described this corre- 
spondence^] on its external analytical side by the transformation [1] 



which is in itself unintelligible, and by equally incomprehensible transition 
from equating to zero a certain expression to the postulation that the space 
integral of the said expression shall be stationary" . 

After more than 40 years after Feynman's statement and more than 80 
years after Schrodinger's admission, there have been many approaches to 
derivation of the law (pQ) from something "we know." Referring only to the 
primary works these approaches have been based on 

a) some ad hoc underlying assumptions, such as the probabilistic nature of 
quantum phenomena [5], or 

b) an earlier analogous assumption about the statistical nature of quantum 
theory [6], or 

1 in Eq.([T} ip{x) is the wave function and H(x',x) is the Hamiltonian of the particle. 
2 namely, between the HJE of a mechanical problem and the "allied" wave equation. 
In Eq. [2] and in what follows, S is the classical action of the particle. 





S = —Lnip 



(2) 
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c) a modification of the classical HJE by introducing into it the dissipative 
term [7], or 

d) the prior knowledge of the Schrodinger equation to be arrived at by mod- 
ifying the classical HJE via introduction of a new postulate (not present 
in classical mechanics, nor following from the experiments) and eventually, 
after a rather involved mathematics, equating - without a physical justi- 
fication - a certain key arbitrary constant to the reduced Planck constant 

m- 

As to the Schrodinger "unintelligible transformation" ([2]) , it has never been 
derived at all. It was used later by R. Feynman as an underlying postulate 
in his groundbreaking path integral formulation of quantum mechanics. 

In sharp contrast to the above works [5]- [8], in this paper 

a) no assumptions about the character of the micro-scale phenomena are 
made, 

b) no new postulates into physics are introduced, 

c) no prior knowledge of the Schrodinger equation is assumed. 

It is nonetheless shown that both fundamentals of quantum mechanics - 
the Schrodinger equation (pQ) and the critical relation ([2]) - are derivable in 
a straightforward (algorithmic) way by using only two well-established and 
well-known rudiments of physics: 

1) The classical relativistic HJE, which follows directly from Newton's sec- 
ond law of motion (see Appendix aEI)- 

2) The experimental fact of energy transfer on micro-scales in discrete quanta. 
2 Massless Particle. 

Let us begin from the relativistic Hamilton-Jacobi equation for a massless 
"particle" [see Appendix A, Eq.(jA13])] 

(f ) 2 " - (3) 

This equation admits two solutions: the particle-like and the wave-like. As 
shown below, these two solutions are not independent but related on a micro- 
scale via the discrete character of energy transfer. This turns out to be the 
"key" which enables us to derive the Schrodinger equation from the above- 
mentioned rudiments. 

3 The numbers of equations in Appendix A have an extra letter A; for instance Eq. 
(All) 
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2.1 Particle-Like Solution. 

For simplicity sake, we consider the 1-D case. A particle-like solution of (1) 
can be obtained by separating the variables: 

S p (x,t) = -Et + f(x) (4) 

where E = constant is the energy of the massless "particle." Upon substi- 
tution of (|3|) in ([3]) and solving it with respect to x, we have 

f = — x + constant (5) 

c 

Without any loss of generality, we can set the above constant = 0. According 
to the definition of the action S, df/dx = dS/dx = p x = p . On the other 
hand, from the energy-momentum relation for a massless particle, it follows 

p = E/c (6) 

Therefore a particle-like solution of ([3j) is 

S p = -Et + px (7) 

2.2 Wave-Like Solution. 

The wave-like solution of ([3D can be arrived at by transforming this relation 
into the wave equation. To this end, we differentiate ([3]) with respect to t 

OS d 2 S _ 2 d^_^S_ 

~dt~d^~ C dx~dtdx~ ~ ( ' 

and then differentiate (|3]) with respect to x 

dS_ j¥S__ c 2 dS_(?S_ _ Q 
dt dtdx ' dx dx 2 

Upon substitution of ([9]) in (jHJ) we get 



d 2 S _ 4 cPS_ Z dS/dx Y 
W ~° dx 2 \ dSjdt) 



(10) 



According to ([3]), we have (dS/dx) 2 /{dS/dt) 2 = 1/c 2 which on substitution 
in (8) yields the wave equation 

^?-4^o pi. 
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For a monochromatic wave, a general solution to (jlip is a traveling wave of 
the frequency oj and the wave number k: 

S w = ae i{ - kx -^\ a = constant (12) 

Thus, by considering one classical equation ([2D, we have obtained two so- 
lutions - relations and (|12p - which describe respectively two patterns of 
behavior of a massless particle: 

• the particle-like behavior, 
and 

• the wave-like behavior. 

As a result, we can say that a massless particle exhibits the so-called 
wave-particle duality feature: on the one hand it behaves like a massive 
particle and on the other hand it behaves like a wave. 



2.3 Derivation of the Schrodinger Relation, Eq.([2]). 

Now we are going to show that the two patterns of behavior of a massless 
particle are not independent. To find the relation between the particle-like, 
Eq. ([5]), and the wave-like, Eq. (|12p . solutions we differentiate the latter 
with respect to t and obtain 

dS M d ( LnS w , 

-uj (13) 



S w dt dt 



On the other hand, 



%=-E (14) 

To relate (fb£|) and (fT3j) . we use Planck's hypothesis and its subsequent 
development by Einstein who has found that the electro-magnetic energy on 
a micro-scale is emitted and received in discrete quanta, that is, 

E = huj (15) 



Thus at micro-scales the energy of a massless particle given by Eq. (j!4[) is 
equal to its expression given by (fT5j) . Therefore by substituting (fT5j) in , 
and using ([13]) we get 

^ = -•, = 1(5**.) (16) 
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Prom (|16p follows 

S p = -Lml: (17) 
i 

where 

^ = (18) 

j4 is an arbitrary constant, and we call if) the wave function. 

Thus we have derived relation (jlT[) directly from classical mechanics (that 
is from "what we know") modified by the existence of the lower bound on the 
amount of energy transfer at micro-scales. This equation is the very same 
"unintelligible" transformation ([2]) "invented" by Schrodinger, and which 
was presumed impossible to derive "from anything we know". 



2.4 Derivation of the Schrodinger Equation, Eq. (fT|) . 

By the definition, the particle's momentum p x = dS/dx . Upon substitution 
in this definition relations (|17p , (JTHJ) , and the wave solution (|12|) we obtain 
the relation between the particle's momentum p and the wave number k 

^=p=*±^ = Wfe (19) 
dx i tp dx 

This is the famous de Broglie relation for a massless particle, which he 
obtained with the help of rather contrived notions of internal and external 
oscillatory processes. If we identically rewrite Eq. (|19p . we obtain the 
relation between the particle's momentum and the wave function in the 
following form 

~^-=W> (20) 
i ox 

Quite similarly, we get from ([14]), CE7|), ([18]) and (|T2j) 



dt il>i dt dt v y 1 

Thus for a massless particle its momentum and energy emerge as the eigen- 
values of the operators {h/i)d/dx and (ih)d/dt acting on the wave function 

As the next step, we differentiate (|20p with respect to x and (|2ip with 
respect to t and obtain respectively 

(hd\(hd\ hdil> a h2 d 2 ^ 2 

\ l OX I \ l OX I I ox ox z 
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d \ /,, d \ ... _ c (t.. B *\ _ rf„,. ... »2#V _ ,.-2 



We rewrite (j22j) and (|23p one more time 

1 d 2 

ip ox z 

= E 2 (24) 



On the other hand, 



£ 2 = (§) 2 (25) 

At the first glance, it seems that two different representations of p 2 and 
E 2 in Eqs. (|24p and (|25p are not consistent. However, this is not the case. 
In fact by using (|17p we have 

_^I^!_, =( ^ ) 2 = _^J_ ( ^ ) 2 

-0 9x 2 <9x ^ 2 ^ ' 

2 1 d 2 dS 2 2 1 2 

Rearranging these equations we obtain 



ipd 2 ip/dx 2 — (dtp/dx) 2 _ d .dip/dx. d 2 ^. ^ 

ip 2 dx ip dx 2 , . 

d 2 4>/dt 2 -(vo 2 a 2 ( ^ 

^ 2 =^W) = ° 

The solution of (EBb 



= - wi) (27) 

gives us the solution (|12p for the wave function of a free massless relativistic 
particle. 
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The derived relations (|22p and f|23|) allow us to obtain the Schrodinger 
equation for a massless relativistic "particle", if we use relation E = cp, 
Eq.©, which follows from the relativistic Hamilton- Jacobi equation ([3]): 

— ft 2 — rip = E 2 xp = c 2 p 2 xp = —c 2 h 2 ——^i/j or 

Ot z OX A /oo\ 

The last equation of (|28p is the Schrodinger equation identical to the wave 
equation (10), multiplied by ft, for a massless particle. 

3 Massive Particle. 

3.1 Derivation of the 3-D Schrodinger Equation. 

Let us consider a free relativistic particle of the rest mass rriQ. The relativistic 
Hamilton -Jacobi equation for a free massive particle ia^| 

{^f-c?{gradS p f = my (29) 

A solution of ([29]) is obtained by the separation of variables 

S p = -Et + p-f (30) 



Inserting (|30p in (|29p we get the special-relativistic energy-momentum for- 
mula 

E 2 = p 2 c 2 + m c 4 (31) 

According to the Lemma 1 in Appendix the substitution identical to 
Eq. (fTT|) for a massless particle 



S p = ALmfj (32) 
transforms equation (|29p into the homogeneous second order nonlinear PDF 



A 2 (^) 2 -c 2 A 2 (gradiP) 2 = m 2 c^ 2 (33) 



4 see Appendix A, Eq.(All) 
5 Equation (B3) 
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where constant A to be determined. Furthermore, according to corollary to 
Lemma 2 equation (|33p is equivalent to the following wave equation: 

A w- Aci d^ + w + ^ ) ~ m ° ciJ=0 (34) 

whose general solution is a traveling monochromatic wave 

ifi = const x exp(ik ■ r — iut) , k(k x ,k y ,k z ), f(x,y,z) (35) 

The respective dispersion relation oj = co{k) is given by Eq.(B15) in Ap- 
pendix B: 

to = ^c 2 k 2 - (m c 2 /A) 2 (36) 

Thus, even a massive particle exhibits the wave-particle duality, formally 
expressed by the existence of two solutions: particle-like (f30l) corresponding 
to the Hamilton-Jacobi equation ()29[) and wave-like (|35p corresponding to 
the wave equation (|3~ij) . 

To find the value of constant A we substitute f)32[) in the definition of the 
particle's energy (fl4|) . use the solution (j35j) . and add the postulate that at 
micro-scales a massive particle's energy is transferred in "chunks" according 
to Eq.CE]), that is 

— E = -frw = = A—-^- = -Aiuj (37 
dt if; dt K ' 

As a result we obtain 

A = - (38) 
i 



which upon substitution in (J32J) yields for a free massive particle exactly the 
same transformation ("invented" by Schrodinger) as we obtained in section 
2, Eq. (fTT|) . for a massless particle. 

Finally, by using (j38j) in the wave equation (f34|) we obtain the Schrodinger 
equation for a free relativistic particle 

1=2 #V ,2 2. 9 V , . . 2 4, n /o Q ^ 

^- nc( ^ + ^ + a? )+m ° c ^ =0 (39) 



From (|38p and (|32|) also follows that relations ()24p derived for a free massless 
particle hold true for a free massive particle. 



^Equation (B21) in Appendix B 
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3.2 Derivation of the De Broglie Relation. 

Now we will derive the de Broglie relation, Eq. (|19|) . for a massive particle. 
By ([30]) and ([32]) we have 

hgradip 

VS p = p = - - — 40 

i ip 

Upon substitution of (|35|) in (j4*0l) we obtain the de Broglie relation 

p = hk (41) 

for a massive relativistic particle, this time rigorously derived from classical 
mechanics supplemented by the requirement of the discrete character of 
energy transfer. This is in contradistinction to the method used by de 
Broglie who employed rather artificial concepts of the internal and external 
oscillations associated with a particle. 

Let us return to the dispersion relation (|36p . By inserting (|36[) in (|35p we 
get 

U) = yj^ + (^-) 2 (42) 

From (|42|) we find the phase and group velocities corresponding to the wave 
properties of a relativistic "particle" at micro-scales 



(43) 



^ = fc =c V 1+( ^ ) >c ' 

du c 2 k _ hk./mo 

~ dk ~ ^ C 2/fc2 + ( moC 2/^}2 = + ( moC 2/^)2 ' 

C 

|v flr | = = < c 

yc 2 k 2 + (moc 2 /h) 2 

On the other hand, the velocity of a relativistic particle is 

p/m /ik/mo 
\/l + {p/moc) 2 a/I + (hk/m c) 2 

where we use (j4"Tj) . Thus in this case the group velocity of the wave coincides 
with the particle velocity of the wave-particle complex. 
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3.3 The Non-Relativistic Schrodinger Equation. 

To obtain the Schrodinger equation for a free non-relativistic particle we 
eliminate the rest-energy from (|39p by representing if) as follows 

4> = Mx,y,z)e- imoc2t/h (44) 
Upon substitution of (|4"4"|) in ([39]) we get 



We divide both sides of (|45|) by c 2 and take the non-relativistic limit c — > oo: 

(46) 



Lim 

c— >oo 



c z or ot \ax z aar oar 



_ , 2 / s 2 a 2 a 2 \ , 

-2thm —- - h —r + —r + — -= )ip = 



which results in 

+ / <9 2 + + _ Q 

at 2mo \ax 2 <9x 2 ax 2 / 

Equation (|47p is the non-relativistic Schrodinger equation for a free particle 
of mass niQ. 



4 Conclusions. 



In this paper the Schrodinger equation for a free particle has been system- 
atically and in a surprisingly elementary fashion derived directly from two 
fundamentals: the classical relativistic Hamilton-Jacobi equation and the 
Planck postulate. 

One of the corollaries of our derivation is the crucial transformation (|17j) , 
lying at the heart of both the Schrodinger approach and the Feynman path 
integral method [10] : 

ip = const x ex.p(iS p /h) (48) 

All the major architects of quantum mechanics Schrodinger [2], Dirac [9], 
and Feynman [TO] introduced (|4"8j) as an ad hoc postulate without any phys- 
ical justification. In our paper this crucial relation emerges as a natural 
consequence of the wave-particle duality. This implies that at micro-scales 
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the wave function tjj replaces the particle's classical action S p , thus reflect- 
ing the dual character of phenomena at such scales. One more important 
corollary is that the canonical representations of particle's momentum and 
energy as the eigenvalues of the partial differential operators are also the 
consequences of the dual character of the phenomena at micro-scales. 

Acknowledgments. The author would like to thank V.Granik for illu- 
minating discussions of the manuscript. 
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5 Appendix A. 

5.1 Derivation of Relativistic HJE from Newton's 2nd Law 

Let us consider the second law of Newton for a relativistic particle in the 
potential field $(r), r(x,y,z) 

dv__J± (A1) 

dt ~ Or [Ai) 
Here the relativistic momentum p of the particle is 

p = — =. v = vv (A2) 

^l-v 2 /c 2 

For the following we express the particle's velocity v in terms of its momen- 
tum p. We obtain from (|A2|) : 

p = p • p (A3) 



a/1 + p 2 /ni^c 2 ' 

Now we consider the momentum p as a function of time t and spatial 
coordinates x,y,z. This means replacement of the Lagrange description of 
motion by the Euler description. In the latter instead of moving with the 
particle, as in the former, a fixed observer watches the different path , which 
the particle can follow at different moments of time. As a result we rewrite 
(|Aip as follows 



dp dp . dp 1 . A 

-77 = jrr + (v • grad)p = -7— + ======= (p • grad)p 

dt dt dt m ^\+p 2 jm\c 2 

H \qrnd.(ri 2 ) — p X cur /pi = —qradQ 

dt 2m v /ITf7^F 

where we use the vector identity 

(p • grad)p = — [grad(p 2 ) - px curlp] 

Since 

^ A , 2/ 2 2^ 1 grad{p 2 ) 

grad(\i 1 + p A jra^c 1 ) = - — — — 

v % m^c 2 y 1 + p 2 / m^c 2 

equation (|A4p becomes 



(A4) 



(A5) 



<9p 2 / ~ TT 1 p x curlp 

— + m c graaU / 1 + p l m^c z - ======= = -grad® (A6) 

dt V 2m yjl +p 2 /mlc 2 
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If we apply operation curl to both sides of (|A6P we get 

d , , s 1 „ px curlp . 

— (curlp) curli „ ) = (A7) 

Equation (1A7I) is identically satisfied, if 

curlp = 

This means that 

p = c/radF (A8) 
where F(r, t) is a scalar function. As the next step we substitute (|A8]) in 



equation (|A4|) and obtain the following equation 



d 

grad(—F + m c 2 \/l + (gradF) 2 /mfc 2 + $) = (A9) 



Its solution is 



F + m c 2 ^l + ( 5 ™(1F) 2 /to 2 c 2 + $ = /(i) (A10) 

where /(i) is an arbitrary function of t. 
We introduce a new function 

S = F(r,t)- J f(t)dt (All) 

Upon substitution of ([ATT]) in (|AT0|) we obtain the relativistic Hamilton- 
Jacobi equation for a free particle of mass mo 

(^ + $) 2 -c 2 ( 5 rad < S) 2 -m 2 1 c 4 = (A12) 
For a free massless particle too = 0, and equation (|A12|) yields 

'dS\ 2 



- 



c 2 {gradSf = (A13) 



dt J 
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6 Appendix B. 

6.1 Auxiliary Mathematical Lemmas. 
Lemma 1 

A nonlinear partial differential equation of the m-th order with n argu- 
ments Xj, j = 1, 2, 3, n 

— dx k j dxjdx k 



dy dy dip 

®'jk...s 

j,k,...,s=l 

is transformed into the homogeneous nonlinear equation of the same order 



^ ay ay oip 



fc=i ft j,fc=i 3 K 



(B2) 



> «ifc...s^l n — a — n — U 

^— ' ax,- axfc ax s 

j,k,...,s=l J 



by the following transformation of the function y: 

y = ALnip (B3) 
where A is a constant and the new function ip = ip(xi,X2, ■■■,x n ) 

Proof 

Let us consider a transformation 

y = y(iP) (B4) 

where ip = X2, x n ). We require this new function y(ip) to be such 

that upon substitution of (|B4p in (|Bip it would transform it into equation 
(lB2l) . Upon substitution of (HU) in ([EE]) we get 



n in, n 



f— ^ dip dx^ jfc eto dxjdxk 

k=l j,k=l 

(B5) 

+ V" a - s (—) m — ^ \ b-0 

^ 3 "' s dip dxjdxk " dx s 

j,k,...,s=l J 
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Assuming that dy/dip / Owe divide both sides of (|B5|) by (dy/dip) r 

El dip ^ 1 dip dip 

, =1 ° fc {dy/diP) m ~ l dx k + ^ fljfc (dy/diP) m ~ 2 d Xj dx k + "" 



(B6) 



E<90 <90 90 & 
jk...S Q /W, fllT. 



dxj dx k " dx s {dy/dip) 1 

Equation (B6) is transformed into equation (B2) if 

1 _ 

or equivalently 

dy_ = A 

dip ip 

where A is a constant. 

Solving differential equation (|B8P we obtained 



(B7) 



(B8) 



y = ALnip (B9) 

The constant of integration has been included into the new function ip. Con- 
versely, upon substitution of (jB9[) in (jBip we obtain equation (|B2|) . 



Lemma 2 

A homogeneous nonlinear partial differential equation of the m — th order 
(to = 2,3, ...) with n arguments Xj(j = 1,2, ...,n) 



E^" 1 ^ + E a jk A^*%L*t + 
fc=i ft j,fc=i 



(BIO) 



90 dip dip 

j,k,...,s=l 

and constant coefficients b, aj,ajk r --,a,jk... s (j,k, ... = 1,2, ...n)has the os- 



+ E ^ m rr-r^ m = 



dilatory solution exp(i ^ aijXj), i = y— T [with the dispersion relation 

i=i 

f(aj,a j: a jk ,...,b) = 0] 
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Proof 

We restrict our proof to a case important in physical applications, when 
the order of homogeneity is m = 2 and the number of arguments n = 4, 
with 3 spatial arguments xix 2 ,xs, and the temporal argument X4 = t. In 
this case Eq. (|B10p written in the Cartesian coordinates is 

a^|i|i + ^ = (Bll) 

CJJb j CJJu 

where the coefficients ajk,b are real- valued and we use the convention of 
summation over the repeated indices. 
Upon substitution of 

t() = const x exp(ia/Xz), 1 = 1,2,3, 4 (B12) 

in equation (|Blip we obtain the dispersion equation /(ay) = where a k (k = 
1,2,3) have the physical meaning of the components k x ,k y ,k z of the wave 
vector k and 04 has the physical meaning of the wave frequency co 

A 2 a^oj 2 + u A 2 [k x (ai4 + a 4 i) + ^(024 + a 42 ) + k z (a u + 043)] + 
+A 2 [(ai2 + a 2 i)k x k y + (013 + a 31 )k x k z + (a 23 + a 32 )kyk z }+ (B13) 
+A 2 (auk 2 + a 2 2k y + 0,33k 2 ) -6 = 

Solution of this equation gives us the dispersion relation to = f (k) for the 
monochromatic wave described by solution (|B12|) of equation (jBlip . 

For the particular case corresponding to the Hamilton-Jacobi equation 
for a free particle we have b = — m^c 4 (c is the speed of light) and 



-c 2 , j = k^4, 

I, j = 4, (B14) 
0, j + k 



The respective solution of (1B13|) is then 

y/c 2 k 2 - (m c 2 /A) 2 (B15) 



Conversely, given function ()B12p and the dispersion relation (|B13P we 
obtain 

dx~dx- k = ^ (B16) 
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As the next step, we multiply (|B16|) by A 2 a,j k , sum the result over all the 
values of j and k ( in our particular case j,k = 1,2, 3, 4) and add the term 

bif; 2 : 

Since by (|B13j) the expression in parentheses of (jB17[) is equal to zero, this 
proves that function (|B3|) is the solution of the original differential equation 
(EH]). 

Corollary. 

Equation (jBllj) is equivalent to the following second order linear partial 
differential equation 

A 2 a jk -^- + biP = (B18) 

CJJb j (J Jb 

with the same constant coefficients as in (|Blip , 
Proof. 

We identically rewrite (|B18p as follows 



2 dip dip 2 _ 2 d 2 ip 
A a jk^—-^ 1- H = ip(A z a jk - — h bip)+ 

OXj OX k OXjOXk 

+A 2 ip 2 a- k \—(— — - — )] = 



(B19) 



ip 2 dxj dx k dxjdx k 
The expression in the last brackets is 

1 .dip dip d 2 ib . d , d r , . d 2 r , _ , 

OXj OX k OXjOX k OXj OX k OXjOX k 

On the other hand, the original equation (|Blip has a solution (|B12p . that 
is ip = const exp(iaixi) . Using this solution in (|B20p we find that 

d 2 

Lnip = 



dxjdx k 



Inserting this back in (|B19p we obtain the wave equation (1B18H . For the 
special case given by (|B14p equation (|Blip is equivalent to the following 
equation 

,2#V A2 2 f° 2 i> , , d2 i>\ 2 4 



A W~ Ac l8# + a# + l#J- m » c> = (B21) 
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